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ABSTRACT 


Th«  flow  of  «  chemically  reacting  gee  through  a  diwergent  noaela 
ifl  ooneldered.  A  pure  diatomic  gee  i«  aesuned  to  diesociate  and  recom¬ 
bine  according  to  eeTen  dletinet  recombination  rate  lawa.  The  effect  of 
coriwotlng  Llghthill's  '•  ideal  dlaaooiatlng  gae"  for  vibrational  oontrl- 
butiona  of  the  molecular  epeciee  and  tranelatlonal  contributions  of  the 
atonic  epeolea  le  shown  to  bo  email.  The  affect  of  recombination  rate 
on  the  flow  of  oxygen^  from  throat  conditions  of  2  atm  and  1^050  K,  le 
shown  to  be  significant  for  reaction  rates  decreasing  with  tenqierature 
and  relatively  unimportant  for  reaction  rates  incrsaslng  with  tempere- 
tuxe.  Bleotronio  analog  computer  solutions  are  presented  for  a  variety 
of  throat  eoayxMitions  and  recombination  rates. 
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A*  PwrloM  Work 

Dmanda  for  laynwlag  ih»  performanoa  of  ehmleal  roekats  naeeaal- 
tata  the  oaa  of  high  anarg^  propaXlanta  with  resulting  higher  tanywratura 
aahauat  produota*  The  analyeia  of  these  high  energy  flow  etroama  baeomaa 
inereaalngly  nora  dlffleult  aa  the  teoparature  la  raiaadi  alnoa  real  gaa 
affaeta  baoooa  nora  prooounoed  and  there  la  aignlfloant  darlatlon  from 
the  elaaaleal  Ideal  gaa  babarlor.  Thermal  phenomena  relating  to  the 
▼Ibrational  and  alaotronie  energy  larala  aa  wall  aa  obamloal  phenomena 
relating  to  the  diaaoelatlon  and  the  lonlaatlon  of  the  different  species 
moat  be  oonaiderad*  The  nta  of  change  of  the  various  flow  paraomtara 
mat  be  erltioally  aaaalnad  to  determine  whathar  or  not  aqulUbrlom  has 
bean  caaehadf  and  avan  moderate  (diangae  auoh  aa  night  be  found  In  the 
expansion  of  a  pure  gaa  through  a  hyparsonlo  nossle  often  result  In  non 
equilibrium  flow  oonditiona.  These  non  equillbrliB  flows  are  particularly’ 
Inportant  for  several  reasons.  In  order  to  examine  the  problem  of  aero- 
dynamic  heating  of  re-entry  vehicles,  the  flow  conditions  of  the  air 
passing  through  tha  bow  shock  nuat  ba  known.  Also,  non  aqullibrium  flows 
auat  ba  ganeratad  in  teat  faoilltlaa  to  almulata  these  re-^ntry  conditions, 
tha  thrust  of  a  roekmt  angina  le  aerloualy  affeoted  Toy  the  extent  to 
idileh  the  propellant  la  dlaaoolated  at  tha  noaale  axlt  and  consequantly 
how  af fleiant  the  engine  hae  bean  In  converting  the  chemicHil  energy  to 
klnetlo  energy. 


Eggsra  (Rafaranoa  1)*  was  oooag  the  firat  to  oorraet  ideal  gaaaa 
to  conpenaata  for  real  gaa  affects.  Ha  aaad  Barthalot»e  aquation  of 
atata  to  account  for  tha  nslaoiaar  else  and  IntaraolaouXar  foreaa  as 
nail  aa  a  Planck  term  to  giw  changes  in  tha  vibrational  heat  oapaoitlaa. 
Ha  than  analjMd  a  nunbar  of  wia  dinaneional  flem  of  a  diatoaio  gaa^ 
finding  that  tbaaa  laparfact  gas  eorraetiona  baoaas  slgoifioant  for  tea- 
paraturas  corresponding  to  tha  flo*  through  a  ahoek  mv  at  a  Mach 
nunbar  of  10.  Chemical  raactlona  mara  not  considered. 

The  problem  of  a  chanicalljr  reacting  gaa  flowing  through  a  nosile 
has  bean  Mamlnad  by  Edae  (Rafaranoe  2).  In  thla  study*  the  gases  mare 
assumed  to  be  in  a  atata  of  perfect  thamodynamle  and  cdianlcal  eqcili* 
brlom  at  all  timas  and  an  effective  iaentroplo  exponent  eas  defined  for 
tha  reacting  mixture,  EquUlhnum  real  gas  effects  have  also  basB  oo»- 
sidared  by  Erickson  and  Creakmore  (Raferenee  3)  and  King  (Raferance  U). 
Ei*lckson  has  examlnsd  tha  flow  of  air  and  pra cants  a  number  of  charts 
for  the  equilibrium  properties.  King  has  studied  the  theoretical  per¬ 
formance  for  an  equilibrium  mixture  of  normal  hydrogen  doling  an  isen- 
troplc  expansion.  Waiter  (Raferenoe  5)  has  included  hsat  sddltion  to  s 
reacting  flow  with  the  hypothaela  that  tha  flow  la  in  thsmcdynsmle 
aquiUbrlum  at  every  point. 

A  Unaarised  treatnsnt  of  reacting  nossle  flow  is  given  by 
Pannsr  (References  6  and  7)*  ^la  aK>roxlmBte  analysis  is  then  applied 
to  tha  two  extrema  caaas  of  vary  fast  reactions*  near  equllllnlum 


*  Raferenoes  are  listed  in  the  Bibliography*  Appendix  A, 
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flow,  and  Tory  Blow  roaotlona,  nar  froien  flow,  Halwa  (Refaranoa  8) 
preeanta  a  dlaouaalon  of  tha  varloua  thaoriaa  eonearoing  tha  osygan 
raeowblnatloD  rata  and  also  an  approxlnsta  aolution  for  ona  diwenaional 
ohannaL  flow,  Datailad  axaot  nuaarloaX  aolutlooa  ara  praaantad  by 
Br^  (Bafaranea  ?),  Ha  conaidared  a  Llghthill  "Idaal  diaaoclaUne 
gaa"  and  aolvad  aavaral  apaclflc  caaaa  aa  wall  aa  preaantlng  an 
approidnata  nathod  for  tha  aolution  of  tha  aat  of  govamlng  aquatlona. 

In  all  axoapt  ona  oaaa  ha  asauaad  tha  raeowbinatlon  rata  waa  aaaantially 
oottstant.  Later  Bray  (Rafaranea  10)  continued  theca  atudiaa  and,  al¬ 
though  ha  dlaoussad  tha  tamparatura  dapandanoa  of  the  reoowblnation 
rata,  oontinuea  to  uaa  a  eonstant  rata  for  noat  of  hia  oalculationa 
(^a  graatly  ainpUfiaa  tha  already  oonplioatad  analyaia).  fraawan 
(Raferanoa  11)  haa  pointed  out  that  ihia  aaauwption  ia  doubtful  but  haa 
praaantad  only  aawi  ooneluaiva  arguaanta.  In  hia  later  paper  Bray  alao 
eonaidara  the  affaet  of  noaala  contour  upon  tha  deviation  from  aqulllbriua 
cowpoaition,  Ba  eoncludad  from  hia  ^iproxljiata  oalculationa  that  there 
are  noaala  eontoora  that  will  have  a  tandaney  to  koap  tha  oonpoaitlon 
oloaa  to  equilibria  but,  in  ganaral,  they  ara  quite  long  and  tharafora 
Impractical,  Additional  analytical  attacka  on  tha  noaala  flow  i^oblaw 
are  given  in  Rafarancaa  12-20  Including  a  nuwbor  of  nuaarlcal  aolutiona 
for  particular  oaaea. 

Wagoner  (Rafarancaa  21-23)  haa  conducted  a  aariaa  of  interaating 
axpariBanta  oonoaming  tha  flow  of  a  reacting  gaa  nlxtura  of  nitrogen 
dlozido  and  nitrogen  tatroxlda  carried  in  nitrogen,  Thia  particular 
siztura  ehangaa  color  aa  tha  reaction  proceeds  (nitrogen  dioxide  ia 
brcwD  and  nitrogen  tatroxlde  ia  colorlcBa)  allowing  concentration 


Maeuresants  to  ba  made  by  light  absorption  tachniqoaa*  Studlaa  of  both 
aaparaonlo  axpanalons  and  flow  through  shook  wavaa  hava  baaa  aadai  and 
tba  raaalta  agraa  nail  with  tiieorotical  predlotions* 

B.  Stataawnt  of  the  Problegi 

This  study  oonoems  the  affect  of  the  rate  of  reooidaixiatlon  on  the 
expanalon  of  a  diatonlo  gaa*  It  viU  bo  aasaioed  that  the  various  onargy 
sodas  are  not  coupled  and  that  the  state  of  the  gas  can  be  daflnod  hy 
the  ideal  gas  law  If  compensation  is  provided  for  the  molaoolar  waight 
changes  as  a  result  of  the  chemical  reaction*  Translational ^  rotatl<»ialf 
and  vibrational  nodes  of  energy  will  be  Ineludedt  however,  all  the  elec¬ 
tronic  energy  levels  will  be  neglected.  An  algebraic  fora  for  the  reoca- 
blnation  rate  will  be  assumed  and  solutions  obtained  for  a  nuaher  of  dif¬ 
ferent  recombination  rata  functions* 


n.  QOTSRHlKa  SQUAT10H3 


A«  On*  TM—ntional  Inyl»cld  Conaenration  Bonatlona 

This  lomttigfttlon  will  consider  cailf  one  dteensionsl  flow  along 
a  atreaiilliie  coordinate  ^  r**  in  a  rl^t  handed^  orthogonal,  curvilinear 
coordinate  ST'atea.  In  the  absence  of  souroas  or  sinks  the  ateady  state 
oontlnnitf  equation  is 

a-0 


f  -  nass  dmsltjr 
A  -  flow  area 
V  -  veloeltgr 

Hegleotlng  eleotronagnetie  effects  and  bodf  forces,  ths  Bonenttn  equation 

la 


dr 


dr 


p  -  presenre 

For  adlabatle  flow  the  enarg7  equation  Is 


4ii 

dr 


^  ^  «  o 
dr 


H  *  enthalpy 


Altbou^  all  sjmbols  are  Included  In  Appendix  A,  Honenelature,  they 
are  defined  as  they  are  introduoed  in  the  text. 
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B.  AiiTtlLarv  Equatipna 
1«  ThenBodynaale  RalatLonahipi 

Only  gM6«  will  be  oonaldared  in  this  stwjyf  and  It  will 

ba  asstaad  that  tbara  la  no  Intaraetionf  or  oottpllng«  batman  tha  vari* 
oaa  anargy  iwdaa.  Tharafora,  tha  partitton  funatlon,  Q,  for  aaeh  dlatoide 
nolacala  nay  ba  arittan  aa  tha  produot  of  tha  indlrldual  partition  fone- 
tlona  for  tha  translational,  rlbrattonal,  rotational,  and  alaetronlo 


BDdaa,  or 


Q  rA  “ 


tv  ft  as 


eV<c 


Sinllarly,  tha  partition  fonetlon  for  aach  atoa  la 

^  ^  traAS  ^^eVcc 

Iha  last  tarn  in  this  axprasaion  sata  the  rafaranoa  lawl  that  la  ehoaan 
for  tha  energy  xaro  for  tha  atonie  apaclaa  at  tha  sam  laml  aa  for  tha 
nolaeular  apadae* 

It  vUl  ba  further  aaaunad  that  tha  componanta  of  tha  partition 
fimetlona  nay  ba  adaquataly  rapraaantad  by 


r  BnMjLll 

Qrot  ’  L  ^ 


QcUc  =  ^ 


l«r*  •*' 


n  -  p«rtiel«  aasB 
Ic  -  Boltsaaim*B  constant 
T  *  aliaolnta  tavparatura 
h  -  Planck' a  eonatant 
Y  -  total  voliana 

^  »  oharaetarlatle  fraquenojr  of  no locular  Tibratlon 
Z  -  BOleoular  nonant  of  Inortla 
(T  *•  ksnmtzy  faotor 

g  -  atatlatloal  wight  (alaetron  doganonej)  of  tho  oloetnmio 
givDad  atato 
E  •>  dlaaodation  onargp 

Tha  partition  funotlona  for  tho  noloouloa  and  atoms  arm  thorofora 


sir'll 


Tho  partition  function  for  a  oolleotlon  of  N  portioloa  is  glwn  by 

q1 


The  factorial  ia  nooeasary  booauao  in  tho  dorlvatlon  of  tha  timnola- 
tional  partition  function  it  was  aasunsd  that  tho  individual  partieloa 
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ware  distinguishable  and  they  are  not.  How,  all  of  the  tharvodynaaie 
functions  may  be  defined  in  tema  of  the  partition  function  and  its  derlwa« 
tlTes.  (It  Is  eonwenient  to  choose  N  equal  to  Avagadxo'a  nua^r  and  use  a 
molar  basis  for  these  calculations. } 

The  HaLofaolta  free  energy  nuqr  be  written  as 

-  kT 

nie  preeeure  and  the  entroj^’  are  related  to  tha  partition  function  by 


The  enthalpy  nay  be  written  as 

H'K+ts  *ipV 

and  tha  Qibbs  free  energy  as 

F  »  /V  -  ^ 

Although  sons  of  the  restrictions  inherent  in  ths  assumptions 
relating  to  the  partition  functions  will  beoons  obvious  or  will  be 
ejEamined  in  greater  detail  in  the  following  sections,  more  extensive  dls* 
cussions  of  these  items  may  be  found  in  References  2l|  and  25 • 

2.  Equation  of  State 

Using  Sterling's  approximation  for  the  factorial  of  large  numbere 


Ih  Nl  ■=  N  V«  M  -  N 

tiM  BaLriioXti  free  margy  may  to  witton 

K=  -RT  [u  Q  .  ,  ] 

and  tha  partial  praaaure  axarted  b7  the  molecular  or  atomic  speelea  la 

•b  --  EL 
V 

Baioritlng  in  tema  of  naae  denalty  and  adding  the  partial  preasurea  of 
the  taD  Blxtore  oaapononte,  the  total  presaura  toeonoa 

In  idtleh  Oc  la  defined  aa  the  naaa  fraction  of  the  dlaaoelated  epeclea 

XnapactlOB  of  Equation  (2^)  ahow  that  It  la  the  olaaaleal  Ideal  gaa 
equation  of  atate  wltb  the  moleoular  Height  corrected  for  the  naaa 
fraction  of  the  dlaaoelated  apadea* 

3*  Bntbalpgr 

Combining  texna,  the  entbalpj  of  either  apeelea  may  be  written 

»  -  «’■'  lir  ^1,  •  "T 

The  anthalplea  of  the  atomic  and  molecular  epeclea  an  therefore 
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•  •!•  R.T  +  ^^E^  ♦  RT 

Sawriting  In  tenu  of  nnaa  onlto  and  adding  tba  contrlbntiona  of  both 
oofoponantai  the  antbalpy  of  tha  Blxtora  la 

Tha  dagraa  of  apprariflatlon  in  tha  ealoalatlon  of  tha  anthalpjr  najr 
be  seen  by  axanlnlng  Plgorsa  2-1  and  2.2*  Uelng  tha  data  of  Halna 
(Raforenea  26} ,  tha  oontrlbatlona  of  tha  aarloof  energy  aodaa  am  aboan 
for  nolaealar  and  atonie  oxygen*  Although  tha  alaetronle  partition 
funotlcm  ban  been  aaaunad  equal  to  Ita  ground  atatlatloal  weight  and 
thua  tha  alaotronio  energy  haa  been  neglaetedf  thla  reatata  In  only  a 
negligible  awount  of  error*  Ihe  electronic  contribution  to  the  enthalpy 
of  Boleoular  oxygen  la  worn  algnlfleant^  aapadally  at  the  higher 
tewparatumai  hoHevari  dlaaoclatlon  deomaaea  the  eonoantratlon  of  tha 
wolaeular  apaelea  at  higher  tewparaturaa  and  oonaaquently  the  alaetronle 
energy  of  the  nolaeulea  nay  alao  be  neglected  without  algnlfleant  exror* 
Alnoat  without  axeeptionf  other  inveatigatom  haye  uaad  the  nodal  of  a 
dlaaoclatlng  gaa  propoaad  by  Lighthill  (Reference  27)* 


60 
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FIGURE  2-1  ENTHALPY  OF  ATOMIC  OXYGEN 
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FIGURE  2-2 


ENTHALPY  OF  MOLECULAR  OXYGEN 
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Equation  (2^)  corrects  Li^thill's  "ideal  dissociating  gas” 
through  the  fact  that  It  alloes  a  variation  In  the  vibrational  ensrg7 
of  the  fflolseular  species  and  the  translational  energy  of  the  dissociated 

species. 


U.  Eqoilibriutt  Constant 


The  eqnllibriua  constant  based  on  partial  pressures^  Tp,  nay  be 
vrltten  in  terns  of  the  Gibbs  free  energy  as 


\v\ 


RT 


(Iha  SDperscript  0  indicates  that  the  free  energies  are  evaluated  at  the 
standard  state  of  one  atnosphere.)  In  terns  of  the  partition  functions, 
the  equillbrlutt  constant  for  the  dissociation  of  a  diatoilo  gas  nay  be 


erltten 


or,  eoabining  the  ideal  gas  law  and  the  partition  functions 

I  -Jm 


I 


b 


-  €<i 


This  equation  represents  the  equllibrlun  constant  only  if  the 
electronic  contributions  to  the  energy  are  negligible.  However,  due 
to  the  high  temperatures  required  for  significant  dissociation,  this 
assumption  is  questionable.  Hanson  (Reference  26)  has  calculated  the 
equillbriun  constant  for  oxygen  including  two  electronic  levels  past 


th0  cround  stat*  for  tha  aolaoula  and  four  alaotroalo  lartla  paat  tha 
ground  atata  for  tha  atom,  Hanaon«a  raaulta  ara  ooaparad  to  tha  valiiaa 
pndLotad  by  tha  abora  aquation  in  Figura  2»3«  Tha  ooaparlaon  la  not 
too  good*  Further  eonaldaratlon  and  an  analjala  of  tha  affaet  of  higher 
alaotronle  atataa  raraala  that  tha  alaetronle  atataa  for  tha  aton  ara 


■uefa  mora  alfftlf leant  than  tha  alaotronle  atataa  for  tha  aolaeula.  Ihla 
fortunata  charaetariatle  alloM  a  moh  battar  ^>pro3dnBtlon  to  the 
aquilibrla  conatatrt  to  ha  nada  by  naglaeUng  tha  rlbrational  tarn  for 
tha  nolaoulaa  to  offaat  tha  alaotronle  imvm,  aaa  Flgora  2-3.  Tharafora 


tha  aqulUbrioa  oonatant  nay  ba  oonaldarad  to  bo 

1/  .  '  1  - 1  ^  ^ 

't  ■  8TV'"'V>  L  N  1  I 


V.T 


In  thla  atudy^  thia  rapwaantatlon  of  tha  aquilibriun  oonatant  haa  boon 
corraotad  allghtly  bar  araluatlng  the  eoafflelant  to  giro  tha  axaet  ra- 
ault  at  a  tanparatura  naar  tha  naximun  tanparaturo  oonaldarad*  Llght- 
hlU’a  "ideal  dlaaocUttng  gae"  uaaa  a  alnllar  approximation  for  the 
aqulUbrlun  oonatant.  Equation  (2-7)  incorporataa  additional  ngor« 
ainoa  It  Inoludoa  an  additional  aquara  root  of  taaporatura  oorraetlon* 


5.  Racombination  Rata 

The  racombinatlon  of  atome  to  fom  diatonic  nolaenlaa  la  ganarally 
eonsidarad  a  temolecular  reaction  In  uhleh  thraa  indlaidual  partlclaa 
participate  in  a  aingla  klnatlo  prooaas,  tha  third  particle  being  naoaa- 
aary  to  renora  tha  energy  of  recombination*  Nunaroua  thaoriea  have  bean 
adaanead  to  predict  tha  value  of  the  recombination  rata  for  thla  type  of 
reaction,  and  Bafarencaa  8,  29$  and  30  contain  dlaoueaiona  of  aotaa  of 
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FIGURE  2-3  EQUILIBRIUM  CONSTANT  FOR  0^  =  20 
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tbaorlat*  Uxafortanatelyt  the  varloiie  epproeohee  do  not 
pertloularlj  regarding  the  teqperetore  dependenoe*  In  genem  It  aegr 
he  that  the  eotieetion  eneror  required  for  the  reocebinetlon 

proeeae  la  negligible  end  the  reoeid)liiAtion  rate  haa  an  algabraio  taa- 
peratnre  dependenee* 

k.  ~  T  "  (.1-6) 

The  Uok  of  asreeMnt  f  roa  thla  point  on  la  evident  qpon  enalnation  of 
lOble  2-1  (abatraeted  froa  Referenoe  29)*  Thla  table  liate  the  velMa 
of  the  reeoidilnation  rate  tar  (Morgan  eorrentljr  aet  forth  In  the  litera- 
tore*  data  for  thla  and  other  high  teaperatnre  reaotiona 

la  aporae  and  qneationable^  ainoa  In  general  the  diaaoclation  rate  ia 
■aaaured  and  the  eqoiUhrlna  oonatant  then  need  to  oalenlate  the  raeoa- 
Unatlon  rate*  The  large  aotivation  energy  required  for  dloaoeiation  ia 
Inelnded  In  an  exponential  tenperature  dependenee  ahicb  affoetivolj  naaka 
the  exact  valoe  of  n  in  Equation  (2-6)*  Thla  atndy  the  affect  of 

varlona  reocnabinatitm  ratea  and  the  foxa  of  Equation  (2-6)  la  aatlafae- 
tozy* 

G.  Reaotion  EaaatA(» 

Diaaooiatlon  of  a  diatonlQ  noleoule  nay  be  oonaidered  to  be  the 
reaolt  of  a  ooUlalon  betueen  a  noleoule  and  another  partiele  if  aoffi- 
elent  energy  ia  avallabla  to  break  the  nDleenlar  bond*  Slnllarlyf 
reoaeblnatlon  nay  be  oonaidered  to  be  the  reault  of  a  eolUaioo  betMon 
two  atone  in  tho  proaenoa  of  a  third  partiele*  The  oppoalng  roeettena 
of  diaaoeiation  and  reoodbinatloD  nay  be  repreaented  aa 
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At  +  B  1  Ak  +5 

ko 

kp  -  dissociation  rats 
•  xsooaMjisticm  rats 
Ths  dlsaacLstion  rats  la  dsf  Insd  as 

i-  —  1  •  tk  f  iiiifill 

MKdtlp  “  L  il  ] 

(Ihs  rats  of  jmduetlon  of  dlsaodUtsd  spaeisa  la  proportional  to  tbs 
eoaoontratloD  of  ths  dlatosiie  apsolaa  tlaaa  tho  total  oanasntratlon  of 
tbs  nlatiirs*}  Tbs  rseoBbinatlon  rata  is  daflnsd  aa 

JL  4Si  I  ,  -  z.k. 

Mk  at  1^  L  I  imk  J 

Ths  not  ehanco  of  ths  atonic  apsolss  is  thsrsforo  ths  bob  of  thoas  tno 
sxpraaalonsa  Bi^andlng  tbs  Biilsrlan  doriratlva  and  rsatrieting  it  to 
stsad^  floVf  this  BOB  bsoonaa 

^  1 
ar  '  V  Mu  I  kfc  i 

k  partioiilar  oass  for  Mbleh  thsra  la  no  nat  rats  of  inrodastleB 
of  sitbsr  obaaloal  apsolaa  la  a  ohanloal  sqnllibrlnii  stats*  For  aqalli* 
brloi  ths  above  aquation  radoosa  to 

ks.  -  il 
ka. 

or^  in  tsma  of  partial  prsaaoraa 


-V 


hz  =  ± 

RT 

Inspection  ahom  that  thla  ratio  of  partial  praaauras  la  aqaivalant  to 
tha  equilibrium  conatant  for  the  roaotlon  being  oonaldared*  Therefore 


kR  RT 


Although  the  above  equation  has  been  derived  for  the  particular  ease 
of  a  system  in  eqallibrlum  (and  thus  no  change  in  ccHnpo8ltlon)|  it  will 
be  used  for  thla  problem  since  it  is  generally  agreed  that  1^,  and 
Kp  are  ell  only  functions  of  temperature  and  not  eompoeitlon  and  the 
relationship  will  hold  for  all  eases  in  idileh  a  mixture  temperature  oan 
be  defined.  Faldnan  (Reference  31)  atatea  that  thla  ratio  nay  be  uaed 
in  non  equilibrium  oases  when  the  individual  partielaB  have  an  eqnili- 
brium  Btatlatical  distribution  of  thair  energy  levels,  nils  therefora 
reduce a  to  the  aaaumption  that  the  relaxation  timaa  for  tho  various 
energy  modes  is  very  short  eompared  to  the  relaxation  tlma  required  to 
attain  chamioal  equilibricn. 

The  species  continuity  equation^  or  reaction  equation,  for  the 
nozzle  flow  problem  may  now  be  written  as 

ar  ■■  V  L  RT  ^  Ma  J 

S.  SuBBarr  of  Ion  DimenslMaliaad  Bouatione 

A  aet  of  eqnationa  deeoriblng  the  flow  of  a  leaeting  dlatcmio 
gaa  have  been  developed.  To  faoilitate  enbeoqoent  oomputetlonB,  the 
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arMf  dlitanoa^  Teloeltjr,  tam»ntar»«  deniitgr^  and  pnoaure  oan  tot  non 
dlmnsionaliaad  vlth  reapaot  to  sona  roformoa  oonditimuif  indioatad  17 
tha  aubaoript  0*  Itaa  raanlting  aquatlona  ana 
Continnitr  Equation 

d  )  =  O  (l-io) 

Monantun  Bouation 


d 


-  ^  dvi 


I 

l  +  «Xo 


Ci-«) 


Sfiargy  Bouation 

vav  ♦DX  V  G(.v-"<n]  =  o  U-'*) 


D= 

R.T. 

(I-IS) 

e=  hi 

V.T„ 

(l-W.) 

U-'ii 

Equation  of  Staty 


T  1,1*  K.^ 

Reaction  Equation 

dx*  ■ 

r'l 

tlMk  4v*^r. 

A  RT,  jC  To'- 

kR 

T  . 

c»  u^ 

\  Kv 

4-^*  T 

l.l-l>) 

ni.  soumoHS  and  discossion  op  results 


A>  General  Con»l«iwatton!l 

referenoe  oonditlone  with  which  the  dleenslonleea  verleblee 
ere  forwed  ere  noet  conwenlentlp  chosen  to  be  thoee  et  the  noesle  throat* 
Exemlnetlon  of  the  eontlnoltp  oQuation  shows  that  this  corresponds  to  the 
point  at  which 

=  -  I 

For  a  wlety  of  particular  cases  this  equation  prorldes  an  expllelt 
definition  for  the  vales  of  ^  (wbioh  Is  a  function  of  the  critical 
vsloelty  at  the  noaele  throat).  However,  in  seneral  this  condition  will 
be  related  to  the  noaals  geanetry  through  the  reaction  equation  and  the 
xeoombinatlon  rate  function*  PhyslcallF,  ^  can  be  considered  to  be  an 
effective  ratio  of  speolflo  beets,  since  It  perfems  a  function  sinilar 

to  that  ratio  in  elesaical  gas  dTnaaies* 

Anelysla  of  the  equations  listed  In  the  preceding  section  shows 
that  both  the  continuity  equation  and  the  energy  equation  nay  be  into* 
grated  directly  with  the  reeulte 


d\J 
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♦  1©  3  I  (VZ) 

Slnco  the  equation  of  state  is  an  algebraic  equation  and  oan  be  solvad 


explicitly,  the  problera  that  remains  is  the  integration  of  the  momantim 
equation  and  the  reaction  equation.  A  careful  study  of  the  reaction 
equation  indicates  that  it  would  be  extremely  difficult  to  Integrate 
analytically,  eren  for  very  simple  eases.  However,  once  the  remainder 
of  the  problem  has  been  solved,  the  reaction  equation  may  be  evaluated 
either  numerically  or  graphically,  to  obtain  the  noasle  contour.  It 
oan  be  partially  integrated  to 


r*  I 


n 


V 


Analytical  solutions  to  the  aomentom  equatlra  ewy  be  obtained  for  a  few 
oases,  but  in  general  numerleal  techniques  mist  be  employed. 


B.  Analytical  Solutions 
1.  .  Proaen  Flow 

If  Uie  reooadiinatlon  rate  is  alow  enough  so  thst  the  compoaltion 
of  the  flow  does  not  change  appreciably  during  the  noasle  expenalon, 
the  flow  is  eaid  to  be  froaen.  This  may  bs  consldsrsd  to  be  the  oeae 
for  wbidi  Bamkohler*a  first  ooeffloient  is 

Oovn  yy  \ 

Sinee  it  is  defined  as  e  chamotenatlo  tlaw  for  reaction  divided  by 


tlw  MBldMiM  tlM  in  tb»  noBsIa.  In  thin  oMt  tte  conditions  at  the 
noaalc  throat  'naat  ba  aueh  that 


_ Z  C  V  -V  ^ _ 

^  ~  (^5  ■*  tt*  \  +  11,  \  ^  ©*■  Jo  Jo  ■*"  O 


t^-4) 


Fignra  3-1  si'roa  this  rssnlt  graphieally  for  tha  initial  eonditioos 
dssoribsd  in  Appandix  C« 

Tbs  gansral  squatlon  for  flow  rslooltgr*  Equation  (3-2)  nay  bs 
spsoialiaad  to  the  eass  of  froasn  flow* 

M^=  +  Jo-  J)| 

Tbs  noaontuB  aquation  ean  bs  intsgratsd  for  this  cast  by  ooiMnlng  it 


with  tha  wslooity  and  ths  sqnation  of  stats*  Ths  rssnlt  is 


2.  Xsothsmal  Flow 

An  isotbsmal  flow  oan  bs  naintainsd  throuf^ut  ths  noasls  sxpan- 
sion  if  snffleisnt  snsrgy  is  islsassd  by  ths  rsoosiblnation  of  atews*  In 
this  oass  tbs  conditions  at  tbs  noasls  throat  wist  bs  such  that 

k^-o  -  - ! -  (?-7) 

+  D  ^ 

This  condition  is  graphically  reprsssntsd  on  Figure  3-1* 

Spscialiaing  ths  vsloeity  to  tha  oass  of  an  isothsnaal  flow,  ths 


rsault  is 
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FIGURE  3-1  THROAT  CONDITIONS 
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i  ,  iilll  - 


Again  tha  Bonentiui  aquation  »ar  ba  Intagratad  if  thie  aquation  ia  oondilziod 
with  tha  aquation  of  atata  and  tha  pvaaauro  nay  bo  aacproaaad  aa 


3,  Equilifariiai  Hoii 

If  tha  reaction  tlm  ia  vary  abort  oonparad  to  tha  raaldenoa  tiaai 


D  o  ^  ^  ^ 


tha  flow  ia  aaid  to  ba  an  equilibrium  flow  and  tha  oompooitlon  ia  a 
unlqua  function  of  tha  aquilibriun  eonatant  and  tha  pvaaaura*  For 
aqnillbriimi  flow  tha  throat  oonditiono  nuot  ba  ouch  that 

d)  *  y-eJO|] 
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fh*  rMtriotion  of  oqiiilifarlaB  flow  dooi  not  prodnoo  «  leant  aiapli- 

floatlon  for  tbo  ■oaantuB  oqoaticmf  and  it  atUl  oannot  ba  Intagratad* 
VoNOrDiia  eharta  and  tablaa  for  oqulllbrliai  propartlaa  hare  boon  poblLshodf 
Roforonoo  3»  to  faeilitato  ntiBBrloal  oaloalatlono. 

U>  lookwrie-Iaoval  Flow 

Althoogh  parhapa  not  too  Intarastlng  In  a  pbfolaal  aanaOf  an 
analjrtleal  aolutlon  nay  ba  obtained  for  tha  oaaa  of  eonatant  praaaara 
flow.  Applieatlon  of  tha  oritarlon  for  tha  noaala  throat  yialda 

,  ivn) 

'  J 

Inapaetlon  of  thla  aquation  ahowa  that  tha  condition  for  a  throat  ra- 
qnlrta  an  asotharnla  dlaaoeiatlon  raaetlon.  If  auoh  a  gaa  wara  iiaad» 
a  itnlqiio  throat  tanparatura  aoold  ba  required  aa  apaelflad  by  tha  abora 
aquation.  A  aomergant-dlTargBnt  noaala  la  lapoaalbla  for  tbo  gaaaa 
eonaldarad  In  Appai^lz  Cj  bouararf  purely  oonrargant  or  dlaergant  aao- 
tlona  nay  prodnoa  eonatant  praaaura  flowa. 

Throuid*  tha  nonantun  aquation  It  oan  ba  aaan  that  a  eonatant 
praaaura  flow  mat  oorraapood  to  a  eonatant  velocity  flow  and  the  nowen- 
tm  equation  radneea  to  a  trivial  identity. 

5.  eonatant  Danaity  Flow 

For  e<Miatant  danalty  flow  tha  etmdltlona  at  tha  noaala  mat  ba 


auoh  that  either 
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or 


Th«  first  rsatriotlon  oorrasponda  to  the  olaeaieal  oaaa  of  ao  laeoapraa- 
aibla  flow  with  no  raaoticn.  The  alternate  restrletlon  reqtULraa  that  the 
throat  compoaltloa  bo  an  eqnllibriaa  oaq>oaltlon«  In  neither  oaaa  la  a 
unique  throat  Talooltf  dafSmd. 


6.  Dleettsslon 

SoM  rather  interesting  results  can  be  dadueed  from  a  oaraful 
atudjr  of  Figure  3-1*  ^  represents  an  effaotlwe  ratio  of  spaelfle 

heats  In  defining  the  orltloal  veloeltF  of  a  reaotlng  gas  flow  at  the 
throat  of  a  noaslSi  For  the  classloal  ease  it  varies  from  7/S  for  a 
dlatonie  gas  to  5/3  for  a  nonatonle  gas*  However^  this  fanlUar  range 
is  Inexaotf  slnos  the  vibrational  energjr  has  not  been  inolnded.  Figure 
3-2  shows  that  ^  inclusion  of  tha  vibrational  tom  lowers  the  effeo- 
tivs  specif  io  heat  ratio  by  approximately  6  per  oent«  Sines  this  la 
the  ratio  of  two  temperature  dependent  funotions  and  varies  less  than 
either  of  then  alone ,  it  could  be  anticipated  that  there  mould  be  a 
considerable  correction  in  the  caleulation  of  the  other  flow  parameters 
if  the  vibrational  energy  is  naglaoted.  The  redeeming  feature  is  that  a 
very  high  prsaeurs  mould  be  required  to  exaotly  approach  this  extrsaa 
ease  at  the  fl»d  temperature  of  li050  X  for  which  the  curves  are  drawn. 
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HoiMTsr^  tha  Inoluaion  of  tha  rlbratioKiBl  anargf  oonttniwa  to  hava  an 
affaotf  of  daeraaalng  Bapiiltadaf  u  tha  eoapooltloa  ohangaa  to  tha  othar 
tattum  of  a  pova  aonatoole  gaa  (oorraapondlng  to  a  aair  low  praoaora). 
Tharafora>  axoapt  for  alaoat  coavlatalj-  dlaaoeUtad  floaa,  tha  vibra¬ 
tional  tama  ahould  ba  oonaidarad,  ngura  3«2  afaoM  that  thia  arror  oan 
axcaad  3  par  oant  in  tha  Taloeity  for  ^  «  0*S0,  Of  oonraa  tha  arror 
will  ba  graatar  for  Baallar  valuaa  of  OC  • 

Tba  for  aqnillbrliiB  flow  natehaa  that  for  froaaa  flow  at  tha 
and  polnta  (tha  Halting  oaaaa  of  tha  aaro  and  Inflnlta  praaaura)  and 
drop#  abarply  in  tha  Intaraadiata  ranga,  approaching  tha  oaaa  of  laothar- 
aal  flow*  Thla  foLlowa  froa  tha  fact  that  In  thia  region  tha  oo^poaltlon 
la  particularly  aanaitiva  to  praaaora*  Tha  valna  of  ^  for  tha  iaothar- 
aal  flow  of  a  raactlng  gaa  la  Indapandant  of  ocaipoaltlon  and  oloaa  to 
uttltgr*  In  tha  Halting  oaaa  of  an  inflnlta  diaaoeiation  anargy*  it  la 
exactly  unity  and  tha  flow  bahavaa  aiallarly  to  tha  olaaaiaal  laothaxaal 
flew,  atudlad  by  Ronar  and  Cadbal  and  reported  in  Rafaranoa  32,  alnca 
no  QiiAwgi*  in  ooapoaltion  would  ba  required  to  aalntaln  the  taBgwratura* 
Tha  throat  raloolty  than  oorraaponda  to  Newton' a  laothanul  apead  of 
sound. 

When  oonaldarlng  tha  flow  of  a  reacting  gaa.  It  la  wall  to  raiwiw 
bar  that  tha  oonpoaition  auat  ba  approaching  aqullllHlaa.  The  deviation 
from  aqulHhrlum  might 'inoraaaa  throu^out  tha  axpanalon  procaas,  but 
thia  nuat  result  from  tito  equilibrium  point  changing  more  rapidly  than 
the  flow  can  follow*  Considering  this,  iaothamal  flow  oan  axlat  wily 
If  (K  >  OC^  t  >lnoa  this  will  naoeaaarily  raault  in  raeonbinatlon. 
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th*  rat*  b*lfig  oontroUad  th*  raooidblnatlon  rat*  fonetionf  F» 

C.  Co*pqt*r  aoitttlon* 

1.  Rcniil*  Contour 

1h*  Ineluilon  of  th*  Bpael**  eontlnulty  •quatlon  with  th*  eonran- 
tlonol  go*  4jrn*Blo  oonramtlon  •quatlon*  ha*  add*d  two  additional  varl«* 
abl0*y  th*  ooqpoaltlon  and  th*  noB*l*  gaoMtiy*  Iharofor*  an  additional 
•qnation  d**orlblng  th*  nosul*  oontour  1*  roqalrad  to  obtain  aolutlon*. 
Th*  prlwaiy  pnziM**  of  this  *tady  1*  to  aarawln*  th*  *ff*ot  of  th*  raeos- 
blnation  rat*  and  tho*  th*  choio*  of  noial*  contour  1*  arbltrarr.  Th* 
■oat  oonranlont  contour 

ha*  b**n  ohoaan  and  1*  shown  In  Figure  >-3,  Tbl*  partloular  ohap*  ha* 
••raral  notewortlqr  eharaotorlatlo*.  PrlwarUy»  it  slnpUfl**  the  equa¬ 
tion*  and  dooraa***  th*  oo^)laxlti**  of  th*  ooqputer  prograa*  In  addi¬ 
tion  It  1*  fluid  Moohanlcally  feaalbl*  sino*  the  nozinni  value  of  the 
half  angle  In  th*  axpanalon  Is  approxlmtely  15  degrees.  In  addition. 

It  foroes  the  reaotlon  aquation  to  predlot  a  raro  darlvatlv*  for  th* 
oonpoaitlon  at  the  throat,  graatly  siBg>IifFing  the  setting  of  th* 
condition*  (see  Appendix  C).  Th*  major  disadvantage  of  this  nosBle  1* 
that  It  doea  not  actually  have  a  throat,  or  nlninon  point,  hut  since  the 
taaperatur*  wist  deorease  as  a  result  of  th*  computer  program,  this 
difficulty  1*  minor. 

2*  Cooputar  Squatlon* 


Considerable  care  must  b*  given  to  the  formulation  of  the  equations 
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WMf  sealing  of  the  variablea  before  thex  ean  be  prognned  for  solution 
on  an  analog  oonpnter  (details  are  presented  in  Appendix  D)«  It  will 
sofflos  to  note  that  sealing  was 

T  •  1  •  lOOr 
p  -  I  •  IOOt 
1  •  lOOv 

t 

j  •  1  •  IOOt 
J  •  1  -  IOOt 

and  the  aaohins  tlas  was  slowsd  down  by  a  factor  of  10.  After  oonsidor- 
able  Bsnipnlationt  the  eoq>uter  equations  ean  bo  written  os  (H  is  the 
reelprooal  of  tbs  area} 
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A  ooaipleta  dlBoossion  of  bonndazy  oondltions  !•  given  in  An>endix 
C.  Three  main  oases  esre  studied  depending  upon  the  initial  oonpositioni 
equUibriuBf  above  equilibrium^  and  belov  equilibrim*  For  eadi  of 
these  main  oases,  seven  subcases  of  different  temperature  dependenoe  of 
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tha  raeonblnatlOD  function  art  aacuaload*  Tha  apaotrun  of  raeoabination 
fonetlona  ara  ahovn  in  Figura  3-Ji. 

3*  Raaulta  and  Dlaoosalon 

aolotlona  fron  tha  analog  oo^mtar  ara  abowi  in  Fignraa 
3^0  and  3*7*  Aa  flow  wariablaa  bahava  aa  axpaetad  and  tba 
praaanro-taaparatura  outra  ia  alightly  boaad  aa  a  raault  of  tha  ineluaion 
of  tba  vibrational  anargjr  aa  wall  aa  tha  raeoadilnatioa  raaotion«  Tha 
valua  of  (aaaantiaUor  tba  of  olaaaieal  gaa  dynanioo)  riaaa  rapidly 
and  than  approaebaa  an  aap^>totio  valua  (for  froaan  flow  thin  noat  lia 
batwaan  li«0  (aonatcaaie  gaa)  and  6,25  (diatonic  gaa))* 

Tba  ehanioal  affacta  in  tha  flow  atraaa  ara  mora  Intaraating  than 
tba  fluid  aaehanieal  variablaa.  Fignraa  3-^,  3-9*  and  3-10  ahow  tha 
change  in  ehaBd.oal  ooapoaitlon  aa  a  funetim  of  temparatura*  Prariona 
inwaatlgatoraj  i»a»|  Rafaranoaa  19  or  20,  hava  indicated  tiiat  at  tha 
low  praaaoraa  atndlad  hara  tha  aaount  of  chaadoal  reaction  oocorrlng 
downatraan  of  the  noiela  throat  ahould  ba  nagligibla*  Thua  tha  aolution 
aaj  ba  adaquataly  rapraaantad  by  tba  froaan  flow  aolution*  fhia  eonoln- 
aion  ia  aubatantlatad  from  thia  atndy  for  tboaa  oaaaa  with  tha  racoidbina- 
tlon  rata  an  inoraaalng  fnnotion  of  tataparatura  (n  >  0).  For  thaaa  oaaaa 
tha  raooBbination  rata  function  bacoiaaa  anall  quickly,  and  thua  the 
ohanieal  reaction  diaa  out  rapidly,  before  a  ai^iif  leant  change  in  oon- 
poaitlon  oan  occur*  Tha  affect  of  the  diaaociatim  energy  la  greatly 
diluted  if  there  ia  only  a  aaall  oonpoaltion  change  over  a  large  tempeza- 
tura  drop*  Since  the  nolaeular  might  correction  ia  a  (1  «  tc  )  tent, 
aaall  changaa  in  eonpoaltlon  will  have  negligible  effecta  on  the 
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FIGURE  3-4 
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NOZZLE  FLOW  PARAMETERS,  CASE  ID 
FIGURE  3-5 
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NOZZLE  FLOW  PARAMETERS,  CASE  HD 
FIGURE  3-6 
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NOZZLE  FLOW  PARAMETERS,  CASE  III  D 
FIGURE  3-7 
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FIGURE  3-10  FLOW  COMPOSITION,  CASE  IK 
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■DlcenlBr  Mlitit* 

Snsprialn^^  dtUbsrat*!/  atartlng  at  tba  throat  with  a  non 
aqailihriOB  ecapoaltloa  has  oitl/  a  anillf  though  prenoonoad,  affaot  on 
tha  oonpoaltlon  profUaa.  In  part  thla  is  a  rasult  of  tha  oholoa  of  tha 
noaala  aontaur*  slnaa^  Ja  this  oasa^  tha  ooapoaitlon  nnat  start  at  tha 
throat  with  a  isro  rata  of  ehanga  in  all  oaaaa*  Thus  tha  nosala  throat 
wmy  oonsldarad  to  bo  at  a  stata  of  psaudo  frosan  onapositlOB. 

Migatiaa  aspononts  (n  <  0)  for  tha  raoodblnation  rata  prodoea  a 
dlffmrant  olasa  of  solutions  for  tha  flow  coapoaitim,  Znspaetion  of 
tha  raarUMi  aquation  will  prorlda  sows  Insight  into  this  dlstinatlon« 
Ibara  an  ***«<*»^^7  two  tans  In  tha  raaotlon  aquation^  ana  indieating 
tha  dirsotion  and  oKtant  to  which  tha  flow  davlataa  from  aquilibrlnn 
and  tha  othar^  oootaining  tha  racawbination  nta  funotion^  indioatlng 
tha  spaad  with  which  tha  raaotion  will  ooour*  For  the  oaaa  of  n  >  0  tha 
rata  of  raaotion  tom  oontrola  tha  idianga  In  oanposltiont  and  ainoa  it 
daoraaaas  as  tha  taaparatura  daoraasaSf  tha  eoaposition  approaohas  a 
oonstanti  or  naarlj  oonstont^  waluo  quieklj.  Ftw  a  nagatiwa  walua  of  n 
tha  nta  of  raaotion  tan  Ineraasas  with  daonosing  tanpantun*  thus 
hawing  a  tandancor  to  inonasa  tha  nta  of  dianga  of  ooqposltion* 

lha  tan  indioatlng  oxtant  of  tha  dowiatlon  fTon  aquilitariun  In 
tha  raaotion  aquation  is  in  ganoral  a  daomaaing  function  with  daemas* 
ing  tanpantun  (although  initially  its  walua  nay  inonasa  dapanding  on 
tha  initial  eonditions)^  and  at  low  tanpantuna  tha  ohanga  in  ooapoai- 
tlon  mot  approach  saro*  Howawar^  for  tha  low  pnasuna  and  pnaaun 
ratios  oonsidarad  han«  tha  tawparaton  doas  not  daoraasa  auffiolantly 


for  this  Moonduy  of  foot  to  jvodoMlmto.  OofortoMtolyf  tho  oo^potor 
OTorloadod  for  t«vioniturM  of  oboat  0»60  for  nogatlvi  a'Ot  and  omm 
ninor  roooallng  of  ^  problm  would  bo  noooooary  to  obtain  oolotloa  for 
graator  tanporatoro  ratloa* 

Tha  oholoo  of  noaala  oontoor  io  of  partiealar  alpilfloanoo  in 
dotandnlng  tha  oonpoaltion.  Iba  llnoar  noiala  otodiad  hara  raotrleto 
tba  alopo  of  tho  oonpoaltion  profila  to  aaro  at  tha  throat*  Ihla 
affootlaoly  aaootha  oat  tha  oonpoaltion  ohanfo  and  dooraaoaa  tha  affoot 
of  raeodbloatlon  rata  on  tha  flow  oarlablaa*  Zt  la  antiolpatad  that 
othar  noaalaa  would  giro  oenaldarablj  dlffaront  oonpoaltion  profUoBt 
partioolarly  for  tha  oaaaa  In  whldt  tho  oonpoaltion  la  not  InltlaU/  In 
aquilibrlon*  Howararf  tho  oonelnoloiia  of  thla  atndf  thonld  not  ba 
affaotad* 

Piguraa  3«12f  and  >13  oonpara  tha  valooltj^  praaaara*  and 

ama  for  tha  two  aatrano  oaaon  of  raoonblnatlon  rata  oonaldarad*  Ao 
wDold  ba  aivaetad,  tha  volooltgr  Inoraaaao  noro  rapidly  In  tha  oonpoal- 
tlon  ohangoB  aa  a  roault  of  tha  diaaooiatlon  anargjr  bolng  rotumad  to 
tho  flow.  Thla  aana  affoot  aooounta  for  tha  faot  that*  for  a  glrwn 
tanparatura  xntlo,  tha  area  la  eonaldarably  dlffarant  for  tha  two  ax- 
trana  oaaaa.  Tha  raault  of  thla  area  Inoraaaa  alao  eauaaa  tha  praaaura 
to  drop  nora  rapidly  for  tha  oaaa  with  groatar  raaetion  (with  reapaot 
to  tanparatura).  Li  tonaa  of  a  glran  nosila  with  a  flxad  area  ratio.  It 
oan  ba  aaan  that  tha  exit  praaaura  la  approxlnataly  the  aona  and  tha 
Ttloolty  la  greater  by  approxlnataly  5  par  oant  for  tho  ooaa  of  nogllglbla 
raaotlon.  Howorar,  thara  la  20  par  cant  nora  thanal  energy  rawaining 
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FIGURE  3-12  FLOW  PRESSURE 
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in  tb»  COM  of  tho  aoro  aotlTO  nijetaro  vhioh  ooaU  potontiAlIj  bo  tvnna* 
lotad  Into  klMtle  onorgr*  In  toms  of  Mob  nnabor^  thin  dlfforonoo  1a 
flow  Toloeltgr  would  bo  oonaidorobly  aoro  pronoimood,  olneo  it  would  In- 
olndo  tho  taoperotoro  oloo* 

All  of  tho  aolutiona  obtolnod  In  thio  otud7  with  tho  oocooptioa 
of  thoM  pMOontod  in  thio  aoetlnof  o*o  Inelodod  In  Appondls  B*  Booulto* 


zv.  oovaoszoiis 


Th»  oonelosiona  of  thlB  Inmatigatlon  0*7  be  eanarlMd  m  follomi 
I.  The  enelog  ooepoter  hu  bean  proven  to  be  e  vereatile  tool 
for  enelirtng  the  flow  of  e  reeottng  gu.  Bo  sejor  chengee 
in  the  progm  need  in  thia  etadf  woold  be  neoeeearjr  to 
erawine  the  eowplete  epeotnai  of  poaeible  boondery  oonditionSs 
2a  Correotions  to  Ligfathlll  * e  "Ideal  dlaeoolating  gaa"  are 
readily  diaeemable  but  not  of  aajor  laportanee# 

3a  The  extent  of  reaction  and  ahape  of  the  ecaipoaltion  profile 
ia  prinarlly  detemlned  by  the  f om  of  the  reaction  rate 

UWa 

Jta  Beaotlon  ratea  that  Ineraaae  with  tanperatuxe  have  little 
effect  on  the  flow  eowpoaitlon  (at  low  praaaore)  and  aueh 
flowa  way  be  aaamed  to  be  froeena 
$a  Reaction  ratea  that  deereaee  with  tewperatura  affect  the 
coagwaltion  algnifieantly  and  in  theae  oaaee  the  ohewioal 
reaction  net  be  eonaidered. 

6a  Initially  non  eqailitarlom  oonpoaitlona  raanlt  in  a  greater 
or  leaaer  obanga  in  eowpoaitlon^  bat  thia  effect  ia  coupled 
very  (doaaly  to  the  particular  noaale  contour  eonaldereda 
7*  The  effect  of  fom  of  the  reootoblnatlon  rate  on  the  flow  of 
a  reacting  gaa  la  algnlflcant  and  waat  be  detemlned  exactly 
before  detailed  analyaea  of  potential  applicationa  can  be  nadca 
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r.  BBCOHMSOMITIOIIS  lOBL  FUTUSS  XMVfiSTKUTOBS 


Tha  faelli^  with  liileh  tha  aat  of  aquatlona  doaoriblng  tfaa  flow 
of  o  xMOtlDg  gaa  Sttggaats  mwit  additional  avaaa  of  Invaatigatlona 

!•  Iha  noBBla  oontonr  baa  a  daflnlta  affaot  on  tha  ooapoaltlon 
profUa  and  a  varla^  of  oontoora  ahoold  bo  oonaidarad*  ftw 
laaa  oonaldorabl7  nora  ooapntlng  squlpnant  than  ws  ooad  In 
this  atadjr  la  amllablaf  thla  aoold  naoaaaitata  ualng  tha 
conaldaftbljr  alapllflad  nodal  of  a  dlasoolatlng  gaa  anggaatad 
b7  LlghthlU. 

2a  Tha  aztanaion  of  tha  atndr  of  a  dlaaodatad  gaa  flow  to  tha 
atudjr  of  althar  a  paralr  lonlnad  flow  or  awan  a  partially 
dlaaoolatad  and  partially  lonlsad  flow  la  raadUy  foraaaabla. 
Tha  graataat  difflonltf  in  thla  aztanalon  would  ba  tha  awaUa* 
blllty  of  flonputar  BOoponanta.  Tha  raaotlon  law  la  naraly  a 
aiopla  loop  and  a  varlaty  of  dlffarant  fora  or  a  aarlaa  of 
dlffarant  raaotlona  oan  ba  addad  aa  aoparata  elroulta  whieh 
can  ba  Inoludad  or  axeludod  at  pradatamlnad  polnta  In  tha 
aolutlox.  In  othar  worda,  aamral  different  nehanlBna  for 
tha  pertlnant  raaotlona,  aaeh  elgnif leant  for  varioaa  proa* 
anraa  and  tanparatoraa,  nay  bo  oonaldarad  atBrnltanaoualy. 

Bray  (RaforHiea  33)  haa  oonaldered  thla  problaa  but  along  tha 
Unaa  of  hla  atudlaa  on  a  diaaocLating  flow,  and  hla  approach 
oan  only  ba  oonaldarad  a  svollninazy  atudy* 
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3«  Ite  nntltlTlty  of  the  oeeposltloa  profile  to  the  reoonbina- 
tien  rote  ■ucgeeta  thet  thle  aer  be  e  naef ol  aethod  for  ezperl- 
aentelly  aeasorlng  the  valoe  of  the  reooabinetlon  rote*  Slnoa 
the  eqnilllrlBa  oonatent  le  well  knoim  end  the  flow  aey  be 
aetabllahid  for  a  raasoiibble  period  of  tlaoi  the  problea  of 
aoMuring  the  dleeooletlaii  mte  and  Inferring  the  reeoaiblna» 
ti<m  rate  te  eliaineted*  Ihe  ooigjooltion  aear  be  aeaenred  verjr 
preoieely  ^etroeoopiaaUj  and  the  low  preaeoree  ooneidered  in 
thle  Bta4r  are  readily  attainable  in  the  laboratory* 
k»  Speoifie  noeele  eontonre  for  particular  applieatione  can  be 
deelpied  ae  an  eatenaion  of  the  oorrant  atudy*  Partieular 
relatlonriilpa  between  the  variooa  paraaetera  oan  be  apeei- 
fied  and  the  reaolting  noeele  eontoor  daterained.  This  analr- 
eia  Bgy  well  indioate  better  epproaohea  for  the  aeaenreaent 
of  moonblnatlon  rate  ttian  that  aoggeated  above* 
in  analog  eoapater  alao  provldea  a  emvenient  aaana  for  ea- 
tbbliahing  dieeontlnoitiea  in  a  flow  field  by  aoddehly  adding 
a  atep  fmetion  (voltage)*  Therefore  it  would  be  oonvenient 
to  aae  it  for  the  atody  of  the  ohenieal  reaotlona  ooeurrlng 
behind  a  atrong  nomal  ahoek  wave* 

6*  The  apead  of  aonnd  in  a  reacting  adatore  ie  a  fonotion  of 
reaction  rate  and  would  lend  Itealf  to  etudy  in  a  iiuinar 
einilBr  to  teat  applied  here*  In  addition,  the  relatlonabip 
between  aonie  veloeity  and  the  oritioal  velocity  at  the 
noBBle  throat  could  bo  exaained* 
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A  -  flow  «rM 

A  -  frM  amrar  (Chaptar  II} 

i2-  typical  dlatoile  wolaoala 

D  -  dlMaaianlaaa  oonatant  (daflnad  by  Sqaatlon  2-15) 

B  -  dlaaoetation  anaiiy 
F  -  (Ubba  fraa  energy  (Chapter  n) 

P  -  dlaenalonleaa  raooeblnatlon  rate  (defined  by  Equation  2-20) 
g  -  atatlatleal  Might  of  the  eleetronlo  ground  state 
0  -  dlMnalmileeB  equUiborltn  oonatant  (defined  by  Equation  2-21) 
b  -  Flanok'e  oonatant 
B  -  enthalpy 

X  -  Mlaonlar  noMot  of  inertia 

i  *  dlMnelonleaa  rlbratlon  ftaoilon  (defined  by  Equation  2-17) 
k  -  BoItaMiin*a  oonatant 
kg  -  dleeoolatlon  rate 
-  reooiri>ination  rate 

Kp-  eqalilbrloB  oonatant  baaed  on  partial  praeeuraa 
m  -  particle  mae 

B  •  ekponant  tor  reombinatlon  rate 

B  -  Arafadro^a  nadMr 

If  -  redproeal  of  the  flow  area 

p  -  preeaare 

Q  •  partition  fnnotion 

r  -  atrMBllne  coordinate 


R  -  onlwrnl  fM  ooiwtant 
3  <•  Mitropgr 

T  -  Absolute  tsnpsrAtors 

V  -  Tsloolt/ 

V  -  total  Toluaa 
Oarsek  Synbels* 

0(  -  dlssooiatlon  panastor  (dsfiosd  ty  Bquatlon  2-5) 

^  -  dlasnaloolsss  constant  (daflnsd  bj  SqaatlcD  2-12) 
-  dlMnslonless  constant  (dsfinad  by  Equation  2-13) 
B  >  dlasnaionlsss  constant  (daf  insd  by  Bquatlon  2-16) 
^  -  eharaetsrlstio  fraquanoy  of  noXsoolar  albratlon 
^  -  naas  density 
IT  -  syuiustry  factor 
Subaorictst ■ 

A  -  atonic  Bpsolas 
D  -  dlBBociatioa 
a  -  aquillbrlun 
slaot  -  alaotronic 

H  -  nolaeular  spseiaa 
R  -  raoonblnaticai 
rot  -  rotational 
trans  -  translational 
Tib  -  Tibrational 
0  -  rsfsrsnos  condition 
SttPsrseriPti 

0  -  standard  stats  of  one  atmosphere 
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ffftbla  C*1  liati  th«  valtwi  of  tbo  ptayiioal  oonitaata  oMd  In  thlt 
Btudj  M  wH  M  Mlooted  proportloB  of  ilx  dlntonie  guoa*  Although 
oxygon  h08  boon  uood  for  tho  epooiflo  oolutlono  prooontod  horo*  tbo  dl- 
nonalonleaB  form  of  tho  oquotlono  oUohb  oono  gonoroliiotlon  of  thooo 
rooulto  to  fit  othor  gaaoo* 

Tho  oholno  of  tho  initlol  toagMrotnro  io  dlffloult  and  dopondo 
opcni  the  oonipoaitlon  and  prooBuro.  More  pertlnant  to  thla  atudy  than 
tho  Initial  tomperatuTO  la  hou  far  tho  Initial  oompoaltion  la  from  an 

oqullltnrluB  congwaltlon*  Using  tho  oubaorlpt  oo  to  Indleato  tho  oquUi* 

% 

hriun  ooapoaltlon  at  the  initial  proaouro  and  tanporatoroB  tho  oqnill* 
hrlon  eonatant  function  nay  bo  written  aa 

^\-  Xoi ) 

Thla  equation  thoa  dofinoB  the  initial  tonporaturOf  through  tho  oqulU- 
brlun  constant,  for  arqr  partloular  oqnlllhrluM  eonpoalti<»  and  Initial 
proaaure« 

Arbitrarily  the  tbIuo  ^  ^  ^  dioaon  to  bo  O.SO.  Tho 

throat  proaaure  has  boon  ohooen  to  bo  2atn.,  thla  low  proaaoro  was 
ohooon  ao  that  tho  solutlona  oould  bo  applied  to  ourront  laboratory 
ayatomo,  nost  of  which  have  glaaa  oalnlng  dianbora  upatroan  of  tho 
noBclo.  For  oxygon,  tho  eorroapondlng  tonporaturo  la  1|0$0K«  Thrao 
dlfforont  Toluoa  of  tho  actual  Initial  oonpoaition  are  studied,  one 
abooe  and  one  below  oquUibriun  and  tho  oquillbrlon  oonpoaition* 

Aa  dlscuBsed  In  Chapter  II,  there  la  oonaldorabla  doubt  aa  to 
tho  oorroct  toIuo  of  the  rocooiblnatlon  rate  for  oa^gon.  Exanlnatlon 
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h  ■  6*62  X  IX)  erg  eec/partlcle 
K  «  6,02  X  10^^  pertlcles/gmHtiole 
R  •  0,31  X  10^  ergs/^K  gmHnole 


k  -  1.38  X  10'^  ergs/^K  particle 

2  2  6  ^ 

1  •  erg  see  /ga  era  •  1,01  x  lo”  erg/cra"’  atm 

TABLE  C-1  Oas  Parametere 
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of  the  Mconbinatlon  rate  function,  P,  shoim  that  the  initial  value  of 
the  recombination  rate  occurs  aa  the  product  of 

is  the  only  place  that  the  characteristic  length  ^)pear8«  Therefore, 
aiqr  change  in  tha  magnitude  of  the  recombination  rate  at  can  be  exactly 
compensated  for  by  redefining  the  value  of  r^.  The  result  la  that  only 
the  temperature  dependence  of  the  recombination  rate  function  need  be  con¬ 
sidered.  The  value  of  F  has  been  chosen  to  be  unity,  a  value  that 

o 

corresponds  to  r^  ■  1.07  <aa  and  kj^  (T.)  ■■  8.U  x  10^  (T^350Cir^  cm^ 
gntHsole*^  800“^  (the  value  predicted  by  Matthens  In  Reference  3U)* 

Seven  integer  values  of  the  temperatura  exponent  in  the  recombination 
rate  were  considered,  -3£n^  3.  Table  C-2  summarlees  the  Initial 

conditions. 

The  critical  value  of  the  velocity,  V^,  ie  difficult  to  calculate. 
However,  for  the  computer  It  is  very  simple  to  determine  its  value  on 
a  trial  and  error  basis,  since  it  appears  only  ttflos  aa  a  ooef  fie  lent 
and,  as  a  result  of  the  particular  noaale  contour  chosen,  only  one  of 
these  coefficients  Is  initially  important,  Thersfore,  ^  la  aet  by 
adjusting  Its  value  until  the  derivative  of  the  area  ie  aero,  the  cri¬ 
terion  for  the  noazle  throat. 
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Tte  aoilvt  altmnt  analog  aoaputarf  or  alaetronlo  dlfforantlal 
analgraari  oonaiata  of  a  eollaetion  of  high  gain  voltage  a^pllflerBf 
reelBton^  eapaoltortt  and  aorvomtera  arranged  in  each  a  oonflgnratlon 


that  the  voltagaa  throughout  the  olronlt  behave  Idantlealljr  aa  the 
phjraloal  parenetera  In  the  differential  equatlona  being  aolved.  The 
elrottlt  la  aat  up  b7  ocaablnlng  a  nuaber  of  Integratere^  aumera^  and 
■ultlpliara  In  a  cloaod  loop  euoh  that  the  operationa  on  tha  voltagaa 
eorvaepond  to  those  required  by  the  differential  aquation* 

The  baalo  elanent  of  the  analog  otapatm  la  a  high  gain  (approzi- 
■ately  l|OOOaOOO)  amplifier.  Thla  aaplifler  haa  tha  oharaoterlstlo  of 
inverting  tha  sign  of  tha  voltage  being  anpUf  led  and  nay  be  aaeoned  to 
have  a  negligible  ourrent  flow  through  it*  Figure  D<01  ahowa  bow  this 
unit  oouplad  with  aerlea  Input  raaistanoaa  and  either  a  feedbaok  reala- 
tanoa  or  eapaoltanoe  ean  be  used  to  add  or  integrate  a  voltage.  Initial 
eonditlona  are  set  by  ia«>OBing  an  InltUl  voltage  on  the  oapaoltor* 
Application  of  Kirohoffa  laws  to  tha  elrouit  of  tha  aoBaar  reaulta  In 
the  output  voltage  being  equal  to 


a 

9, 


«-  -U'-  '-iT.'r.ll 


If  the  aapllfler  galny  A^  la  large  enough  this  nay  be  approzlBated  by 


A  alailar  analysla  of  the  Integrater  and  a  alnllar  approxiaatlon  ahowa 
thatf  in  this  eaeSf 


Ro 


INTEGRATOR 


BASIC  COMPUTER  COMPONENTS 
FIGURE  D-l 
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Thtt  Boat  aignlXloant  faature  of  theoe  tiro  raaulta  la  that  the  ampllfiar 
charactariatiea  do  not  appear  and  the  gain  of  the  circuit  la  a  function 
only  of  the  external  eleotrloal  componenta. 

Another  healo  cnoponant  uaed  in  analog  cooputere  la  the  earvo- 
Bultiplier*  Ihla  unit  eonaiata  of  a  number  of  olroalar  potentionetara 
drlran  fay  a  aervcnotora  One  of  the  potentioaatara  la  oonnaoted  to  a 
referenoe  aivply  and  feeda  ita  output  back  to  tha  aeiromotor*  The 
aermotor  poaltlona  Itaelf  ao  that  the  input  roltaga  la  Juat  balanced 
by  thla  feedback  voltage*  In  other  mrda,  each  follow  up  potentiometer 
la  poaitioned  ao  that  Ita  output  la  decreaaed  in  proportion  to  the  ratio 
of  tha  aarranotor  driving  voltage  to  the  referenoe  voltage.  Thuay  one 
function  may  be  multiplied  by  aa  many  other  functlona  aa  there  are 
follow  up  potentloMtara*  itiltlpllaatlon  by  a  oonatant  ooeffiolent  la 
acooapllahed  by  feeding  the  voltage  through  a  potentiometer  and  tapping 
off  tha  raquirad  paroantaga*  Thla  conatant  eoafflolant  muat  tharafora 
alaeya  ba  laaa  than  unity. 

figure  D-2  ahowa  tha  omoputer  aymbola  uaad  In  thla  report. 

Figure  D»3  ahowa  a  number  of  baalc  oonflgurmtlona  uaad  to  parfom  a 
number  of  oomnon  oparatione.  Solving  for  a^  In  tha  dlvlalon  eireult 


If  tha  ampllfiar  gain  la  aofflolantly  hlghj  tha  aaoond  tern  can  ba 
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DIFFERENTIATION 


BASIC 
FIGURE  D 
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€-  =  -  \OG  — 

**  €1. 

Slmllar3/«  for  ttM  oquar*  root  circuit 
Co  -  VOO 

lbs  diffcrootlation  circuit  Is  only  an  approximation  and,  ao  auch^  is 
gemrally  not  rocofnaendad.  Solving  for  tha  value  of  a^* 

e  =  — ‘  V  (I- a  1 

dt  dt 

For  valuaa  of  a  approaching  unity*  thia  aapraaaion  approxlmataa  tba 
darivativa  mora  cloaaly* 

Tha  alaotronie  analog  conpotar  may  ba  uaad  to  aolva  a  variety  of 
anginaaring  problama  (aaa  Bafarancaa  35-30 )•  It  auitad  for 

aolvlng  ayatoM  of  ordinary  diffarantlal  aqoatlona  daaoribing  initial 
eonditioR  problamaf  hoaavar*  it  nay  ba  uaad  on  non  initial  oonditlon 
problama  or  avan  aoma  partial  diffarantlal  aquationa* 

k  number  of  faetora  muat  ba  carefully  conaldarad  in  preparing  a 
problam  for  aolutlon  on  a  oongxitar*  Aa  avidanead  by  tha  pravioua  die- 
ouaaiMi*  tha  indapandant  variable  for  the  compatar  la  time*  and  all  of 
tba  dapan<tent  variablaa  ara  voltagaa.  Thia  naoaaaltataa  aeallng  of  all 
of  tba  phyaioal  variablaa  In  tanna  of  voltagaa.  Moat  ooiqmtar  oompo- 
nanta  ara  daalgnad  to  oparata  linearly  for  voltagaa  batwaan  0  and  ♦  100 
volta.  Therafore  all  tha  variablaa  ahould  be  acalad  ao  that  thair 
abaoluta  value  la  laaa  than  100  volta  throughout  tha  entire  aolution. 
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(Tho  PACS  eorapatera  oaed  in  this  study  have  an  autonatle  ovarload  indi¬ 
cator  to  atop  the  solution  if  at  any  tine  any  voltage  oxceeds  100  volts.) 
Xn  order  to  alnijiiiee  errors.  It  is  desirable  to  keep  tbs  voltages  as 
large  as  possible,  particularly  in  aubtrmotion  operations.  Scaling  the 
tine  for  tbs  solution  presents  different  considerations.  Long  tines 
allow  errors  to  aoounulate  as  a  result  of  amplifier  drift,  oapaoitor 
leakage,  etc.,  whereas  short  tinea  produes  errors,  sinoe  the  capacitors 
have  their  inherent  tine  constants  and  -Uie  servomultlpliers  will  track 
satisfactorily  up  to  approximately  1  ops.  (For  faster  solutions  non- 
nechanioal  electronic  nultipliers  nay  be  used,  but  they  are  leas  depen¬ 
dable  and  less  accurate  than  servonaltipllere.} 

The  foni  of  the  equations  and  the  sequsnoe  of  operation  is  of 
oonsiderable  inportanos  in  their  solution.  Ideally,  the  equatlone  should 
be  cast  in  a  self-conpensatlng  fom,  that  is,  for  a  decreasing  function 
the  derivative  should  be  equal  to  the  negative  of  the  variable  timss  a 
positive  function.  This  fomulation  has  a  tendency  to  correct  itself 
in  that  If  tha  variable  ia  too  large,  it  correete  itself  by  adjusting 
the  derivative  to  oompeneate.  Effectively  this  ej^roach  foroea  all 
errors  to  oscillate  around  the  proper  solution  rather  than  continually 
increase  aa  the  solution  progresees.  A  serious  difficulty  occurs  if 
the  problem  involvee  the  difference  of  two  values  of  nearly  equal  magni¬ 
tude.  In  addition  to  the  fact  that  mail  errors  are  magnified  in  thia 
operation,  the  character  of  the  solution  nay  change  radically  if  ths 
srror  ia  sufficient  to  tdiange  the  elgebraio  sign  of  the  result.  Thia 
latter  difficulty  may  be  renoved  if  the  nature  of  the  problem  is  such 
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that  tba  corraot  aipi  for  thla  amall  voltage  la  knowia  In  thla  oaaa 
the  result  may  be  forced  to  be  correct  by  feeding  It  through  a  diode 
or  by  driving  a  aervomultipller  and  alloaing  it  to  respond  in  one  direc¬ 
tion  only.  Inapeotlon  of  the  diviaion  droult  shows  that  the  aocuraoy 
is  seriously  affected  if  the  divisor  is  a  very  small  number,  even  though 
the  quotient  may  be  lass  than  100  volts. 

Tba  set  of  eight  simultaneous  differential  equations  describing 
this  problem  of  a  reacting  noazle  flow  have  been  set  up  for  the  computer 
solution  and  are  presented  in  Chapter  HI.  Inspection  of  these  equa¬ 
tions  shows  that  they  are  self  compensating  as  much  as  possible.  The 
reaction  equation  contains  the  difference  between  two  nearly  equal 
quantities,  the  difference  between  them  being  positive  or  negative  de¬ 
pending  upon  the  deviation  from  equilibrium,  and  also,  the  area  function 
starts  at  aero  and  then  Inorsasss.  The  final  progren  is  shown  in 
Figures  C-UA  and  D-liB.  Inspection  shows  how  these  problera  were  circum¬ 
vented  In  this  case.  The  area  function,  lA  dN/dt,  drives  a  ssrvomul- 
tiplier,  allowing  only  negative  values,  sines  the  solution  starts  at 
the  noaale  throat  where  this  function  is  identically  aero.  The  differ¬ 
ence  between  the  two  nearly  equal  quantitlsa  has  not  been  adjusted, 
other  than  aaxlmiaing  the  voltagee,  primarily  as  a  result  of  laok  of 
computer  componante.  The  choice  of  the  reciprocal  of  the  srea,  H,  as  a 
computer  variable  was  made  to  formulate  a  self  ooe^penaating  equation  ae 
well  ae  giving  it  a  maximuD  value  at  the  beginning  of  the  solution. 

Ths  choice  of  pA  «  vsrlabXe  was  mads  to  rsduea  the  rats  of  change 
of  pressure  and  give  more  accurate  reeolts.  The  etaoioe  of  dT/dt  ■ 
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for  the  Independent  verlAble  function  was  made  to  nlnlmlee  the  niuAier 
of  computer  oomponenta  neoeeaaiy  for  the  problem  solution.  Heeerthe- 

I 

I 

leeif  thie  problem  required  the  following  oomponente 

6  -  Integratore 

2li  -  Sueoilng  Anplif  iere 

8  -  Servonultipliere 

(IB-eeparate  nultiplier  eiroulta) 

27  -  Fotentlonetere 

2  -  Switohee 

Figure  ie  a  reproduction  of  a  eemple  eolution  to  the  eet  of 
differential  equatione.  It  ie  the  ou^nit  of  a  aia  channel  Saiftiom 
Recorder*  Despite  the  fact  that  the  program  has  been  optimised  to 
ninSaise  error  and  the  actual  oonponente^  resistors  and  oapaoitors^  are 
acourmte  to  0*1  per  cent  and  eheloead  in  a  themcetatically  oontrolled 
orsn«  errors  will  soouMulato  and  must  be  oheoked.  Althon^  a  ot^lete 
error  analysis  is  out  of  the  question  for  a  problem  this  oomplexi  a 
eery  good  oheck  on  the  accuracy  is  afforded  throuidi  the  fact  that  the 
equation  of  state  may  be  integrated.  In  other  words,  the  error  la  a 
direct  function  of  how  close  the  following  Identity  is  satisfied 

1.  -Av._  .  I 

Figure  D-6  shows  how  closely  the  solutions  obtained  in  this  study  satisfy 
this  restrlotion.  The  curve  presented  represents  the  average  of  a  number 
of  solutions.  Ccaqwrlson  with  Figure  D-5  shows  that  the  aoouraoy  is 
well  within  the  limits  with  which  the  computer  output  can  be  read.  The 
sharp  increase  in  error  at  the  lower  temperatures  is  caused  by  the  snail 
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valuaa  of  noat  of  tha  paraaatera  and  tha  raaultlng  inaoouraoy  of  reading 
the  data.  Hoaarer,  sinoa  tho  njor  intaraat  la  oantarad  upon  the  InltSaX 
portion  of  tha  aolntlon,  thla  preaonta  no  najor  problaa. 
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FIGURE  E-l 
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FIGURE  E-2 
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FIGURE  E-3 
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NOZZLE  FLOW  PARAMETERS,  CASE  IE 


FIGURE  E-4 
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FIGURE  E-5 
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FIGURE  E-7 


86  - 


PRESSURE 

NOZZLE  FLOW  PARAMETERS,  CASEHB 

FIGURE  E-8 
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NOZZLE  FLOW  PARAMETERS,  CASE  JLC 
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FIGURE  E-IO 
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NOZZLE  FLOW  PARAMETERS,  CASE  HF 
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FIGURE  E-13 
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FIGURE  E-14 
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NOZZLE  FLOW  PARAMETERS,  CASE  IHC 
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FIGURE  E-16 
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FIGURE  E-17 


96  - 


1.0  0.5  0.2  0.1  0.05 

PRESSURE 

NOZZLE  FLOW  PARAMETERS,  CASE  MG 


FIGURE  E’I8 


-  97  - 


